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Abstract 

In high energy nuclear collisions, the conventional Glauber model is commonly used to evaluate 
the contribution to J/ip suppression originating from the inelastic interaction with colorless bound 
nucleons. This requires an effective value for the J/^-nucleon absorption cross section which is 
larger than theoretically expected. On the other hand, multiple nucleon-nucleon collisions mediated 
by color exchange interactions, excite their color degrees of freedom. We investigate the importance 
of this effect and find that these excited states provide a larger cross section for J/tjj absorption. 
We conclude that the related corrections are important to explain the effective value extrapolated 
from experiment. 

1 Introduction 

The origin of the anomalous behavior of the J ftp production cross section, as measured in Pb+Pb 
collisions at the CERN-SPS is still debated and several competing interpretations have so far been 
proposed. For a review on the subject see [|], 0]. 

In the early stage of a nucleus-nucleus (AB) collision, the production mechanism of a charmonium 
is controlled by four length scales. They are the coherence length l c < 2 Ej/^/m^^, the formation 
length If ~ 2Eji^ ) /{rn^ } i — m j/^): the mean interparticle distance A ~ 2 fm in a nucleus and the mean 
nuclear size Ra,b — 5 — 6 fm for large nuclei. Their interplay is crucial for the correct theoretical 
formulation of the problem and understanding the nature of early stage suppression. In the kinematic 
regime of the NA38/50 experiment at the CERN-SPS, which detects charmonia with pi a i = 50 GeV, 
one has l c < 2 fm and If ~ 5 fm. Therefore, since l c is shorter than the mean nucleon separation, 
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one can assume that the colorless cc wave packet is produced in elementary nucleon-nucleon (NN) 
interactions. On the other hand, the formation length is of the order of the nuclear size and the 
charmonium takes a finite path in the nucleus before being fully formed. The effect of these length 
scales on hadronic processes in a nuclear environment has been extensively studied ||, |], [|, || and 
experimentally observed in electroproduction of p mesons off nuclei 0. Nevertheless, J/tp suppression 
is mostly treated within the Glauber model, with a constant effective cross section for absorption on 
bound, non-interacting nucleons. Only recently |5|, M, the effect caused by the formation length was 
taken into account in the present context. Notice that while studying the SPS data it is justified to 
neglect the contribution of the coherence length, for RHIC data it will not be so, being If S> l c 3> Ra,b, 
and a new regime will open up. For the present discussion on SPS data, the relevant scales for J/ip 
absorption are // and Ra,b, governing the suppression at the early stage. 

After the early stage, governed by NN collisions, an excited medium with high energy density is 
left behind. It is not yet clear how to describe it, but several points of view have been put forward. It 
is possible that the properties of the excited medium are similar to those of a plasma of de-confined 
quarks and gluons. This was indeed the original motivation for making J /ip suppression so interesting 
0. It has been argued that such a plasma is very opaque to charmonium, therefore preventing its 



final observation as a bound state flO|, 11 1. In the last stage of the reaction, part of the large number 



of produced hadrons, the so called co-movers, is possibly also responsible for suppression |12] , 13, |14}| . 

The present work intends to study the early stage nuclear absorption of J/ip. It is common practice 
to approach the problem assigning a constant value to the J/V'-nucleon absorption cross section. This 
quantity is then considered as an effective parameter, fitted in order to reproduce the production data 
in proton- nucleus (pA) collisions. Our goal is to improve the understanding of the physical phenomena 
behind early stage absorption, neglecting for the time being other possible effects and therefore not 
necessarily being able to describe of the full set of available data. 

Following a known interpretation of the dynamics of AW collisions, we ascribe the main contribu- 
tion to the inelastic cross section to color exchange processes. This assumption is well supported by 
high-energy hadron collision data, showing a plateau in the rapidity distribution of produced particles 
and indicating that the real part of the scattering amplitude is small with respect to the imaginary one. 
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Via color exchange, the 3-quark systems constituting the colliding nucleons become color octet states. 
It was already mentioned in that each colliding system is colored, although the possible effects 
of this were neglected for simplicity. Color exchange may also be accompanied by gluon radiation. 



While the contribution of radiated gluons to J/ij) absorption was treated in |0j, 16], here we focus 
on the first part. Having realized that color-singlet nucleons are found in color octet states after an 
elementary collision, one can expect that repeated scatterings lead to multiple color exchanges. One 
therefore has the possibility to excite the color degrees of freedom of the 3-quark system also to the 
decuplet state. Higher color representations can be obtained when radiated gluons are also present, 
but here we restrict our study on the color dynamics of the constituent quarks of a nucleon. 

In the next section we calculate the cross sections for the interaction of charmonium with a colored 
nucleon. We find an increase of 35 % when the 3-quark system is in octet state, with respect to the 
case when it is a color singlet, while with the decuplet the increase is of 70%. We analyze and discuss 
in detail pA collisions in the third section. Finally, we consider the AB case: the effect of multiple 
scattering is studied and then applied to calculate integrated production cross sections. We find an 
increased absorption, although the experimental suppression in the Pb+Pb cross section cannot be 
explained by the present treatment. The concluding section critically examines the results and gives 
an outlook on possible further improvements. Throughout the paper we indicate with the effective 
charmonium state, composed by a mixture of J/ip, ip' and x- 

2 Interaction of Charmonium with Colored Nucleons 

The main feature we want to stress is the property of the color exchange interaction for the high 
energy collision amplitude. Repeated collisions allow a nucleon to be in any color SU(3) representation 

3 <8> 3 ® 3 = 1 © 8,4 © 85 © 10. To calculate the value of the inelastic cross section for f ona nucleon 
which has undergone multiple scattering, we use the two-gluon exchange model of Low and Nussinov 
17, 18]. Although oversimplified, the model incorporates important properties of high-energy hadronic 



cross sections such as (approximate) energy independence and scaling of the cross sections according 
to the mean square radius of hadrons. One can write at once the absorption cross section, making 
use of the optical theorem, relating it to the forward elastic amplitude. This is convenient, since the 
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Figure 1: Top: Diagram illustrating the meson-nucleon elastic amplitude. The dashed line indicates 
the unitarity cut. Bottom: Vertex T between a gluon and the three-quark system of a nucleon. 

latter amplitude, in the two gluon exchange approximation, is purely imaginary. One obtains 

^ b N = J d 2 r (^,N n \ lY^Vib + Xt-yj)} 2 ^,^). (1) 

i,j 

It is a diagonal matrix element between a product state of charmonium \& and a nucleon N n in 
its color representation n. The variable b is the impact parameter between the scattering hadrons. 
The transverse coordinates 2/1,2,3 refer the quarks in the nucleon, while x\ 2 are those of ^ . In the 
square brackets one has the transverse part of the one-gluon exchange potential. We define its Fourier 
transform as 

v( * - f d ^ 1_ X a exp(ifcz) 1 

where g s is the quark-gluon QCD coupling and fi is an effective gluon mass, introduced to incorporate 
confinement. The index in parenthesis labeling the Gell-Mann matrices indicates to which quark the 
vertex is attached. Light-cone wave functions are used for ^ and N n in the expectation value (|l|). For 
the nucleon, we assume that the transverse positions of the quarks in the nucleon are frozen during 
the early stage of the collision, and therefore use the same spatial distribution both for singlet and 
colored states. 

Inserting in eq. (Q) the expression for the gluon potential given in eq. (|2|) and performing the impact 
parameter integration, one arrives at the relation 

<4ff n = J d2 k ^2^2)2 (tt|l-exp(ifca;i 2 )|fr) 

x (N n \T a (k)T^(k)\N n ), (3) 
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illustrated in the top part of FIG. |T|. The first factor in the integrand is the squared gluon propagator, 
the second contains the phase shift due to the momentum transfer k and is responsible, by the way, 
for color transparency, while the last factor is the interesting part of the present calculation. It is the 
product of two vertices between three quarks and a gluon, averaged with the wave function of the 
colored nucleon. The vertex, illustrated in the bottom part of FIG. |], is 



F a (k) 



16 
3 



+ 4 (A( 2) )^4 exp(iky 2 ) 

+ 44 Wfc' e Miky 3 ) 



(4) 



In order to average the product of the two vertices, one needs to specify the structure of the nucleon 
wave function. The color part of the wave function is the relevant one, while the spatial part is chosen 
to be equal for all color multiplets. We (anti)symmetrize quarks 1 and 2 in order to construct the 
mixed-symmetry octets, therefore use 

l 



^ijk 
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1 

V6 
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J i u j u k 

Because of the SU(3) structure, the different color multiplets have different number of indices. For 
simplicity, we label with n the different multiplets. The factor given by the angular brackets in eq. (|3|) 
can now be evaluated. It involves the contraction of all the indices, leading to 



[N n \T a T^\N n ) = (N\J2 l-c n exp(iky l3 )\N). 



(6) 



The effect of having different color states is all contained in the coefficients c n , listed in the second 
column of TABLE |]. It turns out that matrix elements with n = 8a., 85 are identical. The values of 
c n found, for n = 8, 10, lead to larger absorption cross sections as compared to the singlet case. To 
see this it is useful to define the two-quark form factor F^(k 2 ) = (N\ exp(ife Xij )\N) for the nucleon. 
Since, within our approximation, the latter does not depend on the colour multiplet, we drop the label 
n. Then we can rewrite eq. (Q) as 

a^ n = (*kn(r)i*>, (7) 



where the function in the brackets is 

er n {r) 



& a 2 

3 s 



d 2 k 



{k 2 + fi 2 ) 2 
x [1 - c n F$(k 2 )} 



[1 — exp( i k f)] 



(8) 



exhibiting a similar structure to that of the usual dipole-nucleon cross section, but with the important 
difference carried by the coefficients c n . Notice that with c n ^ 1 and \i = the integral in eq. (||) is 
infrared divergent. Therefore, the cut-off fi, of the order of the inverse confinement radius, is important. 
With fjL ^ it holds that a n (r) — ► for r — > 0, so that one still preserves color transparency, even with 
a colored nucleon. This is indeed a property of the singlet meson. In our work we are not primarily 
interested in the absolute values of cr^- but rather on its dependence on the color state n with respect 
to the singlet value. Yet we estimated the absolute values choosing a s = 0.6 and fi = 140 MeV. Within 
a harmonic oscillator model of the nucleon, one finds that F^(k 2 ) = Ff^^k 2 ), the electromagnetic 
form factor. We assume that this property is general and take the dipole form 

where the last approximate equality allows to fix the parameter A 2 = 12/R 2 , with R p = 0.8 fm. For 
convenience we rescale A 2 = 3A 2 , so that A = 2/R p . With the above definitions, the generalized 
dipole-nucleon cross section can be computed analytically in terms of modified Bessel functions Kq 
and K\, with the result 



cr n (r) 



-DUO [log(0 + Ko(fir) - K (Xr)]} , 
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being = c n /(l - £ 2 ) 2 , D§(£) = 4c n /(l - £ 2 ) 3 , £ = fi/X and H(r) = 1 - zKx(r). The results 

of the calculation, illustrated in FIG. shows the modified dipole-nucleon cross section for the three 



Table 1: Results of the calculation of the meson-colored nucleon cross sections. 
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Figure 2: Dipole cross section for different color states as a function of the quark separation r. Box: 
The quantity A(r) which governs the dependence on color. 

color states of the nucleon. It largely increases when the dipole scatters off a colored object. One can 
get more insight in the effect of the different color states rewriting eq. (||) in the form 



(7„(r) = (r(r)[l + (l-c n )A(r)] , 
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where a(r) is the usual dipole cross section on a color singlet nucleon and the function A(r) is found 
to vary appreciably with r. Its functional dependence is given in the small box in Fig. [2|. To calculate 
the inelastic cross section according to eq. ([?]), we must average a n with the meson wave function. 
Therefore, the dependence of A on r shows that the effect of having colored nucleons is stronger 
for larger size mesons. Recent calculations [19], based on realistic phenomenology for the dipole 
cross section and charmonium wave function, gave values for cross sections of various charmonia on 
a nucleon. At center of mass energy ^/s = 10 GeV, it was found that oj/^tv = 3.6 mb, cr^'N = 12.2 
mb and a xN = 9.1 mb. With a composition of f of 52-60% due to J ftp, 8-10% due to ip' and 32-40% 
due to x> a weighted value of 5.8 mb was obtained. In the present case we evaluate a^jy with a 
Gaussian wave function having root mean square transverse separation (r 2 ) 1 / 2 = 0.25 fm, in order to 
account for the different states (J/ip, ip', x) implicitly included in ^ and therefore adjust the singlet 
cross section to the quoted 5.8 mb. The values found are given in the third column of TABLE [l|. In 
order to minimize the model dependence, we factored out from the colored cross sections the amount 



corresponding to the singlet case, therefore writing cfj^y- = 5 n . The values for 5 n are given in the 
last column of TABLE H. They constitute the important result of this calculation , showing that \& 
scatters with an color octet nucleon with a 35 % larger cross section with respect to a singlet nucleon, 
while with the decuplet the increase is even of 70%. The large values found now require a careful 
analysis to establish the validity of this newly proposed absorption mechanism. 

3 Analysis of proton-nucleus data: is the color mechanism already 
at work ? 

Since the absorption pattern observed in pA collisions is usually considered to be the baseline for the 
AB case, it is necessary to examine it in some detail in the light of the possibility of a color mechanism 
for charmonium absorption. 

Consider the illustration given in FIG ^. There one can see the longitudinal and transverse 
geometry of the collision between a proton and the row of target nucleons with which it collides. On 
the left hand side one sees the projectile P which collides with the target nucleon Ti, producing a would 
be ^> meson. Successively the projectile interacts with the target nucleon T2 which, then, collides with 
the VP meson. According to our initial argument, each nucleon-nucleon collision leaves the nucleons in 
a color octet state. As illustrated on the right hand picture, since ^ is produced by P, it is overlapping 
with it. If then the charmonium is struck by T2, P and T2 must also overlap, indicating that they have 
previously interacted. Therefore, in a pA collision all the nucleons with which ^ has an interaction are 
color octets. This means that the effective absorption cross section normally used in Glauber model 
calculations in order to reproduce pA data, should be the cross section with a color octet nucleon. On 
the other hand, another simple geometrical argument shows that this might not be the full truth. In 
fact, in order to claim that VP interacts with a colored nucleon, the charmonium state itself is required 
to be well separated from the parent projectile 3-quark system, preceeding ^ in propagation through 
the nucleus. This takes a finite time, since the charmonium is produced inside the projectile and has a 
smaller velocity. In the rest frame of the nucleus, the difference in velocities between and the parent 
3-quark system is Av ~ m\/2E^ — m 2 /2E 2 , being v = yl — m' 2 /E' 2 ~ 1 — m 2 /2E 2 . Therefore, one 
requires a finite length L = R p / Av, with R p = Rp /7 P , in order to have the charmonium separated 
from the projectile. In the following we focus on the kinematics of the NA38/50 experiments at the 
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Figure 3: Left: Longitudinal space-time representation of charmonium production in a pA collision, 
in the cm. frame of a NN collision. Right: Illustration of the overlap between projectile (P), parent 
target (Ti) and nucleon with which scatters inelastically (T2). 



CERN-SPS, for charmonia with laboratory momentum pq, = 50 GeV and consider projectile protons 
with at least E p = 200 GeV. Then, since 7 P = E p /m p , taking Rp^ = 0.8 fm we obtain L = 2 fm. The 
physical meaning of L is the following: After its creation, the charmonium has to travel a length L 
(in the target rest frame) before one can safely assume that the next interaction occurs with a colored 
nucleon. Since the mean length which \& travels inside a nucleus is 2 /3 Ra , the retardation length L 
reduces the effect of colored nucleons considerably. It destroys it in light nuclei with A < 50, and it 
reduces it to about 50% in heavy nuclei. 

Now, before we analyze the data for charmonium production in pA collisions, we discuss the 
importance of charmonium formation time, introduced in a quantum-mechanical way in |,§. Within 
a simplified version of a two-channel model ||], the ^N absorption cross section evolves in the eigentime 
r of ^ as 

ofSr(r) = 4W + (^ (0) " 4V) cosCr/r/) , (12) 

where tj = (m^t — mj^)~ l = 0.3 fm is the charmonium formation time, <7^ N is the inelastic 
cross section and a (°) = af^{T = 0) is the cross section of the so called "pre-meson" on a nucleon. 
As mentioned in the last section, calculations give <t'^ n = 5.8 mb Unfortunately has not yet 
been calculated. We therefore use the value 2.7 mb from Ipfl. 



Making use of eq. (12), one can then write the expression for <!' production in pA collisions, nor- 
malized to one nucleon, as 



pA 



NN 
0> 



d 2 bdz 



(13) 
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We now want to compare this expression with the conventional one obtained in the Glauber model, 



containing only an effective absorption cross section, therefore making the replacement cr|^(r) — * a\ 



In this way one can perform the z-integration and obtain the well known formula 

. 1 



JpA ^NN I Mi 

<J\<, = cr-s, / a b 



l-a e ^T A (b)/A A 



rr eff A 



(14) 



being T A (b) the nuclear thickness function. 

For A = 208 and py = 50 GeV, using a realistic profile for Pb |2(| and forcing the two expressions 
m eqs. © and (@) to be equal, we obtain cr^y = 3.8 mb. This number has to be interpreted as 
the value of the effective absorption cross section of on a color singlet nucleon, including formation 
time effects. To obtain the corresponding value for the case when ^ scatters on a color octet nucleon, 
according to the result of the previous section we have to enhance the singlet result by 35%, obtaining 
5.1 mb. The values of = 3.8 mb and 5.1 mb are indicated as full lines in FIG. |||. Because of a 
possible retardation of the color mechanism, we expect the experimental values to lie between these 
two lines. 

We have then re-analyzed the existing data taken by 3 fixed target experiments: NA3 at CERN 
PHI , which measured the ratio of VP production cross sections between Pt and H2 for 200 GeV projectile 
protons; NA38/51 |2^| at CERN, which measured the cross section for different targets (H2, D2, C, 
Al, Cu, W) struck by 450 GeV protons; E866 at FERMILAB ||], which measured the ratio of 
cross sections between W and Be for 800 GeV protons. In order remove the effects on absorption 
due to different ^ momenta, we impose the same kinematic conditions, considering VP's of 50 GeV 
momentum in the laboratory frame, corresponding to xf = 0.15, 0.01, —0.04 respectively, for the 
three experiments. We then assume that the absorption mechanism is parameterized by the effective 
cross section <r^y, appearing in the Glauber model expression given in eq. (0). The ratio of cross 
sections is defined as 

R A /B = <Ty A / ■ (15) 

Therefore, while the pA cross section depends on the two parameters a^ N and c^y, the ratio depends 
only on the latter. We made three separate fits using the MINUIT package from CERNLIB, including 
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quoted statistical and systematic errors for the various cross sections and ratios. We extract the 
effective absorption cross section by minimizing \ 2 ■ The result of these fits is presented in the second 
column of TABLE ^. The average is obtained, using the weights Wi = l/Saf, asa = J2i w i a i/J2i w ii 
while the error is 8a = (J2i w i)~ 1 ^ 2 ■ The values extracted from experiment are represented as open 
points in FIG. |j. The mean value of 3.5 ± 0.6 mb seems to favor the singlet mechanism, for which a 
cross section of 3.8 mb is theoretically expected. 

On the other hand, the assumption that is produced on a bound nucleon with the same cross 
section as on a free one, used in eq. (|lj|), can be an oversimplification of the physics under study. 
We know that in a nuclear environment parton distributions are substantially modified [^4], |25|, pfj. 
In particular, the gluon distribution in nuclei exhibits anti-shadowing effects at XBj — 0.1, i.e. the 
gluon distribution gA(x2,Q 2 ) in a nucleus is larger than the one in a nucleon g(x2,Q 2 ). We expect 
that charmonium production on a nucleus is enhanced by a factor R 9 A (x2, Q 2 ) = gA{x2, Q 2 )/g{x2, Q 2 )- 
We estimate these effects by rescaling the \I / cross section with the A-dependent correction R 9 A = 



0.032 log(yl) + 1.006, extracted from 2C]. For a heavy nucleus this amounts to an enhancement of 
the production cross section by a factor of nearly 1.2, therefore requiring a larger cr^jy. We fit the 
data once more, obtaining the results listed in the third column of TABLE |2|. One sees that the new 
average value of 5.7 ±0.6 mb now seems to favor the octet mechanism, for which we estimated a value 
of 5.1 mb. The results of the fits, together with the theoretical expectations, are all plotted in FIG. |], 
where the full points are the result of accounting for anti-shadowing. One observes that the situation 
is, in all, quite uncertain. 

The above analysis is based on the assumption that there is no dependence of nuclear effects 



Table 2: Results of the fits to pA data for the various experiments considered. See text for details. 



Experiments 
p f = 50 GeV 


*Tn [ mb ] 


Glauber fit 


Gluon corr. 


NA3 


4.5 ± 1.7 


7.3 ± 1.9 


NA38/51 


7.1 ± 1.6 


10 ± 1.8 


E866 


2.7 ± 0.7 


4.7 ± 0.7 


Averages 


3.5 ± 0.6 


5.7 ± 0.6 


Theory 


3.8 {1} 5.1 {8} 
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Figure 4: Results of the fits to pA data for the various experiments considered. Open dots refer to the 
Glauber model fit while full dots account for gluon anti-shadowing corrections. The horizontal lines 
indicate the theoretical expectations for the singlet and octet absorption mechanisms, respectively. 

on the proton energy, is the charmonium energy is the same in all experiments. This is justified in 
many phenomena, e.g. the coherence time effect Q. However, there might still be some dependence 
on the proton energy. According to p7|], energy loss of the projectile gluon in the nuclear medium 
substantially reduces the production rate of charmonia at large xf and is still rather important at 
small xf. In order to compensate for energy loss, the gluon distribution must be evaluated at a higher 
value of xi, therefore suppressing ^ production. This will in turn reduce the final state absorption 
on nucleons in a fit to the data. To account for this effect, already observed in Drell-Yan data |28| ]. 
further study, which goes beyond the scope of this paper, is required. 

To summarize, the data seem to suggest that ^ is absorbed by colored nucleons in pA collisions. 
The picture is though not yet conclusive due to the large experimental errors. Nevertheless, we follow 
our theoretical prediction and examine the potential role of the color mechanism in the more complex 
case of AB collisions. 

4 Color dynamics in nucleus-nucleus collisions 

Preparing the ground for the forthcoming discussion, it is useful to recall the conventional expression 
within the Glauber model for the production cross section at impact parameter b. Generalizing 
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eq. (13), to the case of AB collisions, one simply obtains 



* - (6) = J d 2 b A dz A d 2 b B dz B 6 2 (b A + b B -b) a§ 



NN 
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A-l 



B-l 



d 2 b 

x p A (z A ,b A ) [l-a%%T A (z A ,b A )/A 

x PB(z B ,b B ) [l-af s N T B {z B ,b B )/B 

x S F . s .i.(z A ,z B ,b A ,b B ) , (16) 

which describes the collision of two rows at impact parameters b A and b B , overlapping due to the 
(5-function. p A and p B are the nuclear densities while 

r z A 

T A (z A ,b A ) = dz' A p A (z' A ,b A ) and 

J — oo 



_, f+OC 

T B (z B ,b B ) = dz' B p B (z' B ,b B ) 

J Zrt 



'zs 

are the thickness functions. The last factor Sf.s.i.( z A, z b, b A , b B ) represent the probability of survival 
due to any additional "anomalous" suppression factor which goes beyond conventional nuclear effects, 
that is final state interactions with the produced matter (QGP, hadron gas, etc.). Here we are not 
primarily interested in these contributions, therefore set Sf.s.i. = 1- On the other hand, we focus on 
the reinterpretation of the terms describing nuclear effects discussing, in the next two subsections, color 
exchange in the context of multiple scattering, then applying the obtained results to the calculation 
of ^ production in AB collisions. 

1 Multiple color exchange in a nucleus 

Having completed the computation of the VP cross sections with colored nucleons, we need to look at 
how these color objects develop in a collision between nuclei. Consider a nucleon of one of the nuclei 
colliding with a row of nucleons of the other nucleus. It undergoes multiple scattering, with repeated 
color exchanges. Keeping the projectile quark coordinates 1/1,2,3 frozen, according to the high-energy 
approximation, it is possible to calculate the probability V n ({yj},z,b) for a nucleon in color state 
n = 1,8,10 after having traveled through the row up to z at impact parameter b. This can be done 
in a very elegant way, by solving an evolution equation for the color density matrix of the 3-quark 
system. The problem has been considered some time ago with the general aim of studying the color 
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exchange interaction in the nuclear environment. A detailed calculation 2£] shows that 



Vn 



27 



1 + 20.F + 2 g {1} 
16-40 F + 8 Q {8} 
10 + 20.F -10 Q {10} 



where 



and 



Z J) = exp [ - ^ ^ er(yy) T(z, &) ] 



(17) 



(18) 



<?({%•}, z,b) = Y: exp[-| (7( Wi ) T(z, 6) ] . (19) 

The coefficients T and are expressed in terms of the standard dipole cross section cr(y^), appearing in 
eq. flll|), which one can construct by pairing two of the three quarks in a nucleon. Their arguments are 
the differences of quark coordinates yij = yj — yj. The thickness function T(z, b) can be taken as that 
of nucleus A or B, depending on which nucleus one is considering. We can average the probabilities 
with spatial nucleon wave functions as 



P n (z,b) = (N\ V n ({yj},z,b) \N) 



(20) 



therefore integrating over internal quark coordinates. Due to the structure of eq. (O), the averaging 
procedure expressed by eq. ( |20| ) is transferred to an average of T and Q. In analogy to the definition of 
P n (z, b), we will define also the functions F(z, b) and G(z, b). This average can in principle be exactly 
evaluated, although in practice one needs to know the correct form of the dipole cross section and of 
the spatial distribution of the three quarks in a nucleon. We therfore choose to average the exponents 
appearing in eqs. (|l~8| ) and (|l9|), instead of the full exponentials, i.e. evaluating exp(— (N\X\N)) instead 
of (N\ exp(— X)\ N) for a given X. This corresponds to neglecting inelastic shadowing corrections 
which, according to |?0], pHl , are small. Due to rotational symmetry one has (a(yij)) = 2/3 a 
all We took the value cr™ N = 30 mb. In this way one obtains 



^ivfor 



1 



1 + 20 F(z,b)+ 6 G(z,b) {1} 

P " ' ~27 X \ 16 " 40 + 24 G ( z ' $ { 8 i 

10 + 20 F(z,b) -30 G(z,b) {10} 



(21) 



where 



F(z, b) = exp 



9 „in 

'8 



(22) 
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Figure 5: Evolution of the color probabilities Pi, P$, Pio, along z and with b = 0. Also shown is the 
profile ppb(z, 6 = 0) of a Pb nucleus. 



One notices several properties of the probabilities. First of all, the limit T — > implies F,G — > 1. 
This means that Pi — ► 1 and Ps,Pio — > 0. In other words the nucleon, before entering the nucleus, 
is in singlet state as it should be. Moreover, Pi + P$ + Pio = 1 for any T, therefore probability 
is always conserved. Finally, if the nucleus is large enough one has T 3> 1 and F,G -C 1. This 
implies that if sufficient scatterings have taken place, the color probabilities reach the statistical limit, 
given by the first coefficients of eqs. (plj). In other words, the 3-quark system becomes unpolarized in 
color after many collisions. Out of 27 states, one is the singlet, 16 are in the two octets and 10 are 
in the decuplet. The ^-dependence of the color probabilities is illustrated in FIG. ||, together with 
the longitudinal profile ppb(z,b) of a Pb nucleus. Soon after the nucleon has penetrated the nuclear 
profile, a process that involves ~ 4 fm corresponding to 2 nucleon mean free paths, the statistical limit 
is essentially reached in which about 16/27 ~ 2/3 of the nucleons are in octet state while 10/27 ~ 1/3 
are in decuplet. The amount of singlet is negligible. 

With the calculated probabilities P n and the previously evaluated color cross sections cr|^ , we 
can now construct the cross section for ^ when it scatters with a nucleon, taking into account all its 



and 



G(z,b)=ex V \-\a i S N T{z,b) 



(23) 
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color states. We obtain 

a%%(z,b) = J2 «$k p n(z,b). (24) 

n 

As obvious from FIG. ^ the projectile nucleon which collides with several target nucleons quickly 
reaches a "color equilibrium" . Using the values for S n given in TABLE |, that is 1.35 for the octet 
cross section and 1.7 for the decuplet one, the effective absorption cross section of f on a nucleon, 
within a nucleus-nucleus collision, becomes 



1 16 10 
— + — x 1.35 + — x 1.7 
27 27 27 



x <J^p N 



b=0 

~ 1.47 x offi = 5.5 mb, (25) 



where the value cr^v = 3.8 mb has been used, according to our calculation in the previous section. 
We can therefore summarize saying that, although the cross section for with a singlet nucleon is 3.8 
mb, including formation time effects, its effective value in AB collisions with large nuclei is 5.5 mb. 

As a final remark, notice that the calculation leading to eqs. (^T|), (p2|), ( p3|) and ( |24| ) is performed 
assuming that each projectile nucleon undergoes collisions with singlet nucleons in the target. This 
neglects the color excitation of the target. A complete treatment of the problem is very involved and 
is not expected to change the picture drawn. In fact, as seen in FIG. ||, color excitation is a rather fast 
process and accounting for the neglected features would make it even faster. Moreover, as discussed 
in the introduction, gluon radiation, here neglected, would speed up color excitation even more. 

2 Charmonium Suppression in AB collisions 

Having discussed in detail colored cross section and multiple color exchange, we obtained, with eq. (^ 
the important ingredient needed to calculate in detail the production cross section in AB collisions. 
It is first of all necessary to discuss some kinematical features concerning ^ absorption and follow 
the longitudinal space-time diagram of FIG. ||, where the kinematics is illustrated. Considering the 
absorption caused by one of the nucleons of nucleus B, with coordinate z'b, one notices that, before 
interacting with the meson at point Db, it has traveled through nucleus A from +oo down to za- 
The position za depends on the velocity in the center of mass frame. This, in turn, depends on 



xp = 0.15, from which one obtains that = Jx F s/{Am\, +x F s), where s is the center of mass 
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Figure 6: Longitudinal space-time representation of a row on row collision taking place at impact 
parameters (b A ,bB) within a AB collision. The production point is indicated with P, while the 
absorption points by Da and Db, do to nucleons of A and B respectively. 

energy of a NN collisions. With an elementary geometrical reasoning, making use of the triangle 
PQDb for nucleus B, and DaCbP for nucleus A, one arrives at the relations 



z A = z A - (z' B - z B ) , z B = z B - (z' A - z A ) R^ 1 



(26) 



where i2\j- = (1 — vy) / (1 + vy) . These must be used as arguments for the effective cross section in 



eq. (24). Therefore one has two cross sections, one for absorption by nucleus A, taking into account 
the multiple scattering of its nucleons with nucleus B, and one for nucleus B, depending on A. We 
can now modify eq. (|l(i| ) by making the replacements 



(27) 
(28) 



af^ T A (z A ,b A ) -» / dz' A Eb(zb,1>b) PA{z' A ,b A ) , 
cr^>N T B (z B ,b B ) -> / dz' B T, A (z A ,b A ) pB(z B ,b B ) ■ 
The position-dependent effective cross sections 

Z A (zA,b A ) = °$N n Pn(~ZA,b A ), (29) 
n 

^B(~z B ,b B ) = £ ag s Nn P*(zB,b B ), (30) 

n 

are obtained from eq. (p4|). 

Now eq. (16) can be used to evaluate total cross sections, with a further integration in impact 
parameter. We compute values to compare with all available data measured by the NA38/50/51 
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collaboration at the SPS |32(| . We use realistic nuclear profiles |20(] and scale the center of mass 
energy according to the parameterization a§j N ~ (1 — mj/^/y^) 12 . The results are compared with 
the data in FIG. |. The dashed curve corresponds to a standard Glauber model calculation with the 
singlet absorption cross section G^jy^ — 3.8 mb, as evaluated in Sect. 3, while the thick curve refers 
to the improved calculation which takes into account color excitations, studied in Sect. 4. In order 
to emphasize the difference between the two absorption mechanisms, we neglect here anti-shadowing 
effects. For the pA calculation we increased the value of cr^^ by 35%, therefore used =5.1 
mb. This is the extreme situation, when ^ always encounters color octet nucleons. For the AB case 



we used eqs. (29) and fl30|). Since our main intent is to compare the the new absorption mechanism 
to the conventional Glauber model one, we do not present the calculation of the -Ex-dependent cross 
section for Pb+Pb collisions, which we cannot explain. 

Different from the Glauber model result, our improved one provides a better agreement with the 
data, except for the Pb+Pb point. In other words, the proposed absorption mechanism due to color 
dynamics is found to be an important contribution to nuclear effects in ^ production. It must be 
accounted for, together with other mechanisms of suppression |l5|, [H], |33| which we here put aside. 
In the search for new physics one must in fact rely on a solid baseline, employing the best of our 
knowledge of the dynamics of ^ absorption at the early stage of a nuclear collision. 

5 Conclusions 

In this work we studied the dynamics underlying ^ suppression due to nuclear effects, employing the 
idea that color exchange processes in multiple NN collisions may play a significant role. 

First of all, we calculated inelastic cross sections for ^ scattering on colored nucleons and found 
an increase of about 35 % when the nucleon is in a color octet state as compared with the singlet case, 
while the enhancement factor for the color decuplet is about 70 %. 

Making use of available calculations of ^fN cross sections, we included formation time effects and 
obtained the effective value for ^ absorption on singlet nucleons crj^y =3.8 mb. Scaling this value 
with a 35% increase we obtained o"|f^- g =5.1 mb and with a 70% increase f|^v 10 = 6.4 mb. 

We then performed a fit to available data on ^ production cross sections in pA collisions. The result 
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Figure 7: Calculations of total J/ip production cross sections in AB collisions and comparison with the 
NA38/50/51 data. The curves, drawn to guide the eye, correspond to Glauber absorption (dashed) 
and to absorption by colored nucleons (full). 

of the fit, including gluon anti-shadowing corrections, is that the effective absorption cross section for 
is cr^y = 5.7±0.6 mb, therefore suggesting that in pA collisions, the produced charmonium interacts 
with color octet nucleons. The conclusion is, though, rather weak because the experimental points are 
quite scattered. 

Total production cross sections for ^ were then calculated and compared in FIG. [7| with available 
data for AB collisions. Compared to the Glauber model result, the agreement is improved and we 
conclude that effects due to color exchange processes are substantial and provide the value = 5.5 
mb for the effective absorption cross section on nucleons in AB collisions (See eqs. ( p4|) and (|25|)). 

In conclusion, we have completed a study of nuclear effects in charmonium production in AB 
collisions, providing a better baseline for future calculations. These amount to an improved treatment 
of gluon bremsstrahlung, including nondinearities arising because of gluon cascades, to the calculation 
of effects caused by gluon energy loss and to the evaluation of Sj_-dependent cross sections. This will 
hopefully provide the starting point from which to address more precisely the nature of anomalous 
suppression in Pb+Pb collisions. 
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